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AN  INTEGRAL  EQUATION  FORMULATION  OF  THE 
EQUATIONS  OF  MOTION  OF  AN  INCOMPRESSIBLE  FLUID 


INTRODUCTION 


In  the  course  of  various  investigations  it  became  apparent  that  it  might  be  pos¬ 
sible  to  express  the  equations  of  motion  of  an  incompressible  fluid  solely  in  terms 
of  integral  equations.  In  fact,  it  turns  out  to  be  possible  to  derive  such  a  set  of  cou¬ 
pled  integral  equations  in  what  may  be  called  the  vorticity-velocity-enthalpy  formu¬ 
lation.  This  report  contains  a  derivation  of  these  equations  and  a  discussion  of 
their  properties. 

The  formulation  employed  here  is  similar  to  that  used  by  Howe  (reference  1) 
in  his  investigations  of  acoustic  wave  equations.  In  that  work  Howe  employed  a 
formulation  in  which  the  acoustic  pressure  was  replaced  by  a  stagnation  enthalpy 
that  included  the  dynamic  pressure  term.  The  present  formulation  generalizes  the 
incompressible  form  of  the  stagnation  enthalpy  integral  equation  to  include  the 
viscous  term  and  formulates  the  rest  of  the  equations  of  motion  in  terms  of  inte¬ 
gral  equations  as  well.  Various  mathematical  manipulations  are  then  carried  out  on 
these  equations  to  render  them  in  their  final  form. 


MATHEMATICAL  FORMULATION 


Consider  the  flow  of  an  incompressible  fluid  in  both  bounded  and  unbounded 
domains.  The  velocity  may  then  be  expressed  as  the  sum 

Uoo  +  U,  (1) 

where  in  an  unbounded  fluid  domain  u  is  the  disturbance  velocity  and  Uw  is  the 
freestream  velocity  (U„  =  U„(t)),  and  in  a  bounded  fluid  domain  u  is  the  total 
velocity  and  U„  may  be  set  to  zero  (see  figures  1  and  2  for  illustrations  of  bounded 
and  unbounded  domains).  The  governing  differential  equations  are  then,  the  conti¬ 
nuity  equation 
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Figure  1.  Bounded  Fluid  Domain 


c 


V 
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and  the  Navier-Stokes  equations 

d(U^+u)  +(uw+u) .  y  u  -  -  +vV2m.  (3) 

The  Navier-Stokes  equations  may  also  be  written  as 

+VB  -(U00+B)Xtt=-vVxtt,  (4) 

dt 

where  the  specific  stagnation  enthalpy  B  is  defined  as 

B-Rl£2  +  i[(0.+n)-(U.+o)-U_-Uj,  (5) 

P  4 

and 

n  (6) 

tt=  V  xu 

is  the  vorticity. 

The  usual  conditions  on  the  flow  boundaries  are  the  no-flux  boundary  condi¬ 
tion 

n  u  =  nuB  on  S,  (7) 

and  the  no-slip  boundary  condition 

n  x  u  =  n  x  uB  on  S,  (8) 

where  S  is  the  boundary  of  the  fluid  domain  V,  and  uB  is  the  velocity  of  the 
boundary.  For  unbounded  flows  there  is  also  a  condition  at  infinity  that  may  be 
expressed  as 
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and  G  is  any  vector  Green’s  function  of  the  form 

G=r+H(r),  (12) 

where  H(r)  is  a  regular  vector  function  and  r  =  x  -  § .  This  identity  holds  for  any 
vector  field  that  is  differentiable  and  for  which  the  integrals  exist.  For  reasons 
that  will  become  clear  subsequently,  it  shall  be  assumed  that  the  curl  of  G  is  always 
zero. 

The  second  integral  identity  is  a  generalization  of  Green’s  third  identity  (refer¬ 
ence  3): 
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where  G  is  any  scalar  Green’s  function  of  the  form 

G  =  i+H(r),  (14) 

and  H(r)  is  a  regular  function.  This  identity  holds  for  any  scalar  field  that  is  differ¬ 
entiable  and  for  which  the  integrals  exist. 

The  integral  equations  will  be  derived  for  the  case  of  a  bounded  flow  domain. 
If  an  unbounded  flow  domain  is  under  consideration,  care  must  be  taken  when 
considering  the  contribution  of  the  flow  boundary  “at  infinity.”  In  that  case  the 
integrals  over  the  “surface  at  infinity”  may  be  eliminated  by  assuming  that  the  dis¬ 
turbance  velocity  and  the  vorticity  vanish  sufficiently  fast  at  large  distances  and  by 
defining  the  stagnation  enthalpy  as  in  equation  (5)  so  that  it,  too,  goes  to  zero  at 
infinity. 

From  the  integral  identity  (10)  and  the  fact  that 


V-«=  V-(Vxu)  =  0, 


(15) 


one  immediately  finds  that 

rr 


rrr 


/3u  =  -  #[<  n  •  u)G  -  G  x  (n  x  u)]dS  - 


jj 

s 


Gx»dV, 


jjj 

v 


(16) 


which  expresses  the  velocity  field  in  terms  of  the  vorticity  in  the  fluid  domain  and  a 
boundary  contribution.  This  equation  is  merely  a  generalization  of  the  Biot-Savart 
law.  Similarly,  equation  (3)  may  be  rewritten  as 


V  x  «  = 


_  vbI 

at  J  ' 


(17) 


Therefore,  the  integral  identity  (10)  yields 

/3te  =  -  (jj)  [G(n  •  to )  -  G  x  (n  x  to  )]dS 


+  1 
V 


rrr 


(18) 


G  x  .  (u„+u)  x  to  +  VB]dV. 
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However,  since 


VB  x  G  =  V  x  (BG)  -  BV  xG, 
=  V  x  (BG), 


(19) 


where  it  has  been  assumed  that  the  curl  of  the  vector  Green’s  function  G  is  zero, 
one  finds  that  the  enthalpy  is  required  only  on  the  boundary  of  the  fluid  domain, 
and  equation  (18)  becomes 


j8»  =  -  [«(  n- 


)  -  Gx(nx«)]dS  +  - 


1  BGxn  dS 


s 

+ 


1  f  Gxf(U^  -  (Uw+n)  x  tt]dV. 

JJJ 


(20) 


Note  that  the  time  derivative  may  be  expressed  as  a  material  time  derivative,  yield¬ 
ing 


rr 


/3®  =  -  (fj>  [G(n-  to)  -  Gx(nx©)]dS  +  i  <H)  BGxndS 

s 


JJ 

s 


if 


'}}  “f3! 


^  -  \  V[(U.+U)  -(U„+u)]}dV;  (21) 


or,  operating  on  the  gradient  term  in  a  manner  similar  to  that  employed  in  equa¬ 
tion  (19),  one  finds 

/3to  =  -  4  [G(n  •  to)  -  Gx(nx  to)]dS 


JJ 

s 


+  I  C)  [b  \  (U„+u)<(U„+u)]Gx ndS 


JJ 

s 


+ 1 
V 


rr 


JJ 

v 


Gx  D(U-+a>  dV. 
Dt 


(22) 


(See  Hildebrand  (reference  4)  for  the  applicable  vector  identities.) 
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Finally,  to  obtain  an  equation  for  the  specific  stagnation  enthalpy,  note  that  if 
the  divergence  of  equation  (4)  is  taken,  one  is  left  with 

V2B  =  V-[(UC0+n)x®].  (23) 


Similarly,  equation  (4)  can  be  employed  to  show  that  on  the  boundary 


9B 

9n 


js  =  n  •  [  -  -(U”+U)  +  (U.  +  u)  x  «  +vV2u], 

=  n.[.a(Uo°+u)  +  (Uoo+tt)Xtt.vVxtt]. 


(24) 


Hence,  employing  integral  identity  (13)  and  using  expressions  (23)  and  (24)  derived 
above,  one  finds  that 


/3B  +  ®BffdS=<tl> 


JJ 


|-n •  +  n •  [(U^-t-u)  x  o ]  -  vn-  (V  x  »)j-GdS 


rr 


V-[(Uw+u)x«]GdV. 


(25) 


Employing  the  identities 


(V |  x  »)G  =  Vg  •  («G)+VGx  tt, 

Vj  -(ux«)G  =  V|  •  [(ux  «)G]+VG-(ux  «), 


(26) 


and  the  theorems  of  Gauss  and  Stokes,  one  can  find  the  third  and  fourth  terms  on 
the  right-hand  side  of  equation  (25)  yield 


•  V  x  <a)GdS  =  -V 


jj 


rr 


dp  n*  V  x (ttG)dS  -  dt)  n*(VGx«)dS 


jj 


=  -v 


k* 


rr 


a>G-dS-dt)n-(VGXtt)dS 


JJ 

s 


rr 


=  v  dp  n  •  (VG  x  )dS, 


jj 


(27) 
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and 


I 


rr 


{v-[(Uro+u)xtt]}  GdV  = 


V-{[(UTO+u)x©]G}dV 


VG-[(Uw+u)x«]dV, 


n-j^U^+i^x  «]G}dS 


rrr 


VG-[(Uw+u)x«]dV. 


Therefore,  equation  (25)  becomes 


jj 


/3B+#B^dS 


at 


•(VGx  to)jdS 


rr 


VG-[(U„+u)x«]dV, 


(28) 


(29) 


or,  since 


n-(VGxtt)  =  -VG-(nx  «), 


(30) 


one  finally  arrives  at  the  integral  equation  for  the  enthalpy: 


/3B  + 


B^dS 

an 


|n*  8^+  --G  +  vVG •  (n x  a ) jdS 


rrr 


VG-[(U„+u)xte]dV. 


(31) 
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Therefore,  the  set  of  coupled  integral  equations  in  the  vorticity-velocity-enthalpy 
formulation  is,  for  the  case  of  general  boundary  motion, 


and 


and 


/3u 


=  -  <j)  [(n •  u)  G  -  G x (n X u)]dS  -  J 


Gx  »dV, 


Cf  Cf 

/3w  =  -  (!)  [G(n-tt)  -  Gx(nx «)]dS  +  U.)  B  (Gxn)dS 
JJ  JJ 


+i 


rr 


G  x  .  (u„+u)  x  tt  jdV, 


rr  cc 

/3B  +  a  )  B  —  dS  =  -(  )  jn- d(U^)-  G+ vVG -(nXM )]dS 


JJ 

s 


JJ 

s 


+ 


rrr 


VG-[(U„+u)x«]dV. 


If  the  free-space  vector  and  scalar  Green’s  functions 


G  = 


(32) 


(33) 


(34) 


(35) 


and 


G  =  - 


(36) 


are  employed,  then  the  integrand  of  the  time-derivative  term  in  equation  (33)  may 
be  written  as 


rx|t(D.+u) 
r5  " 


=  v(i)x|t(U.+u), 


(37) 


so  that 
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f/rrxfi(u»+'*) 


r  rri 

Jj 


dV  = 


rrnx-(V„+u)  rrr 

f—r - dS  jjjf(F)dV- 


(38) 


Therefore,  with  the  Green’s  functions  given  by  equations  (35)  and  (36),  the  set  of 
coupled  integral  equations  (32)-(34)  may  be  written  as 

/3u  =  -  J>  [fn'")r  -  liiiHiliOlclS 
r  w  1  r3  1*3 


i 

V 


rx  to 

r3 


dV, 


(39) 


and 


rr  r 


/?$&  —  “  dp 


(n-  te)r  _  rx(nXtt) 


r3 


rr 


dS  +  idb  B 


(rxn) 


r3 


dS 


and 


+  i 

V 


T  nX-|(UM+u) 


dS 


s 

rrr 


1 

v 


J  (1)  +  r  x[(Uqo+u)  x  »])  dV 
\  at  V  r3  J  ’ 


(40) 


/3B+ 


B  |.(I)  dS  =  -  (jp  /n  •  2£k±H)  ffl  +  v  r  S*  x 1}  dS 
5nV  \  at  r3  / 


r'[(Ua+u)  x  te]  HV 
r3 


V 


(41) 


In  two  dimensions,  the  free-field  vector  and  scalar  Green’s  functions  are 


G  =  £  , 
r2 


(42) 
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and 


G 


>n(i), 


(43) 


respectively.  The  corresponding  integral  equations  in  two  dimensions  are  then 


=  -  rx(nxu)jH1_ 

s  JJ 


(44) 


and 


*•-  f  +\\\  h  I^K) 


dS 


I  ff  J  3*  In  (I) 
v  [  3t  W 


rx  [(U^H-u)  x  to 


r2 


ljdS, 


(45) 


and 


+  f  B I ln  ©  dl  =  ■  f  {“ '  I (D" +u)  ln  (f) + v  £t1}  dl 


+  f  r-[(Uw+u)Xtt]  ^ 

JJ 


(46) 


where,  in  two  dimensions, 


2 tt  in  V 
ir  on  S. 
I  0  in  Vc 


(47) 
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DISCUSSION 


The  boundary  conditions  expressed  in  equations  (7)  and  (8)  may  be  readily 
implemented  in  equation  (39)  by  substitution  in  the  integral  over  the  boundary, 
yielding 


/3u  =  - 


f(n  •  ttB)r  _  r  x  (n  x  uB) 


dS- 


rxa 

r3 


dV. 


(48) 


Equation  (40)  also  includes  an  integral  over  the  boundary  that  involves  the  dot  and 
cross  products  of  the  surface  normal  and  the  surface  vorticity.  The  treatment  of  the 
boundary  terms  here  is  not  as  straightforward.  It  is  shown  in  appendix  C  that,  for 
rigid  body  motion  with  a  no-slip  condition,  the  dot  product  of  the  normal  and  the 
vorticity  can  be  expressed  in  terms  of  the  boundary  condition;  that  is, 


n  •  to  =  -2n  •  ft, 


(49) 


where  the  boundary  condition  has  been  decomposed  into  its  translational  U  and 
rotational  O  parts  as 


ub  =  U+r  x  ft. 


(50) 


However,  the  cross  product  of  the  normal  and  the  vorticity  cannot  be  expressed  in 
terms  of  the  boundary  conditions  since  it  involves  derivatives  of  the  velocity  field 
normal  to  the  boundary  (see  appendix  C).  Thus,  for  rigid  body  motion,  equation 
(40)  becomes 


/3to  — 


if 


r(-2n  •  ft)  _  r  x  (n  x  to ) 


r3 


d  S  +  I  <0>  B 


(r  Xn) 
r3~ 


dS 


+ 


jj 

s 


nx|(U„+u) 
- - - dS 


rrr 


v 


rx[(Uoo+u)x  to]\HV 

r3  J 


(51) 
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Note  that  equation  (41)  involves  the  same  cross  product  boundary  term  and  there¬ 
fore  cannot  be  simplified  by  using  the  boundary  conditions. 


Equations  (48),  (51),  and  (41)  now  represent  the  system  of  equations  to  be 
solved  both  in  the  interior  of  the  fluid  and  on  its  boundary.  However,  the  nature  of 
the  equations  changes  depending  on  whether  or  not  the  boundary  or  interior  is 
being  considered.  In  the  interior  of  the  fluid  the  equations  take  the  form 


47711  =  -  (H) 


JJ 

s 


(n  •  uB)r  _  r  x  (n  x  uB) 
r3  "  r3 


rrr 


dS- 


r  x  « 

r3 


dV, 


(52) 


and 


rr  r 


4tt  a  =  -  dp 


r(-2n-Q)  rx(nXtt) 
r3  ’  r3 


dS  +  i 1 


(rxn) 


dS 

r3 


+  i(0) 


dS 


JJ 

s 

rrr 


1 (1)  +  rx[(U,+u)x»] |  d Vj 


(53) 


where  each  equation  has  been  written  with  its  unknowns  on  the  left-hand  side. 
Note  that  no  equation  is  needed  for  the  stagnation  enthalpy  in  the  interior  of  the 
fluid  since  equation  (53)  requires  the  enthalpy  only  on  the  boundary. 

If  one  writes  the  equations  for  the  boundary  values  in  the  same  form  one  finds 


rr 


rr 


27TOJI  -  dp  **i“*»l  dS  =  -  qp  r(-2n'°)  dS  +  ljbB^dS 


JJ 

s 


rr 


+  -  <n> 

v  w 


nx#t(u“+u) 


jj 

s 


dS 


s 

'  tfffffi  (f)  +  rX[(-*-rr)-<j>--}dV'  (54) 
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and 


rr 


+ 


jjj 

v 


r-[(U.+tt)x»]HV 

r3 


(55) 


Note  that  no  equation  is  required  for  the  velocities  on  the  boundary  since  they  are 
already  specified  by  the  boundary  conditions. 


Equations  (52)  through  (55)  represent  an  integral  equation  reformulation  of 
the  equations  of  motion  of  an  incompressible  fluid.  These  equations  possess  some 
interesting  properties.  First,  they  contain  no  spatial  derivatives.  Second,  they  only 
require  knowledge  of  the  “pressure”  quantity  (the  specific  stagnation  enthalpy)  on 
the  boundary  of  the  fluid  domain.  Third,  if  the  variables  u,  a,  and  B  are  consid¬ 
ered  as  independent,  then  these  equations  are  linear  in  each  variable,  a  condition 
one  might  call  “pseudo-linear.”  Fourth,  since  all  the  volume  integrals  present  in  the 
equations  contain  the  vorticity  and  since  the  farfield  boundary  condition  is  incorpo¬ 
rated  in  the  formulation,  it  is  immediately  apparent  that  the  domain  of  interest  in 
the  fluid  may  be  restricted  solely  to  that  region  of  the  fluid  in  which  the  vorticity  is 
nonzero.  These  facts  naturally  lead  to  speculation  as  to  the  usefulness  that  this  for¬ 
mulation  might  have  from  a  computational  point  of  view.  The  fact  that  the  equa¬ 
tions  contain  no  spatial  derivatives  suggests  that  they  might  not  be  prone  to  the 
requirement  of  artificial  viscosity,  which  plagues  so  many  finite-difference-based 
computational  schemes.  Similarly,  since  the  integral  equation  formulation  only  re¬ 
quires  knowledge  of  the  “pressure”  on  the  boundary,  it  may  obviate  the  need  for 
elaborate  schemes  to  guarantee  the  convergence  of  the  pressure  calculation  in  the 
interior  as  is  required  in  many  differential  approaches.  The  “pseudo-linearity”  of 
the  equations  may  allow  the  use  of  efficient  iteration  schemes  for  the  solution  of 
the  equations.  Finally,  the  ability  to  restrict  the  domain  of  consideration  to  only  the 
domain  of  nonzero  vorticity  should,  in  typical  high  Reynolds  number  cases  where 
the  vorticity  is  exponentially  small  outside  of  a  thin  region,  greatly  reduce  the  num¬ 
ber  of  nodes  at  which  the  solution  must  be  obtained. 
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SUMMARY  AND  CONCLUSIONS 


An  integral  formulation  of  the  equations  of  motion  of  an  incompressible  fluid 
has  been  derived.  The  equations  consist  of  a  generalization  of  the  Biot-Savart  law 
for  determining  the  velocity,  an  integral  expression  of  the  momentum  equation  for 
determining  the  vorticity,  and  a  boundary  integral  equation  for  determining  the 
stagnation  enthalpy.  The  equations  are  linear  in  each  independent  variable,  with 
the  nonlinearities  entering  only  through  cross  terms  of  the  vorticity  and  velocity. 

This  formulation  possesses  several  salient  features,  including  the  total  absence 
of  spatial  derivatives,  the  fact  that  the  stagnation  enthalpy,  or  pressure,  is  required 
only  on  the  boundary  of  the  fluid  domain  and  the  fact  that,  since  the  vorticity  is 
present  in  all  volume  integrals,  the  domain  of  integration  in  this  case  is  restricted 
to  the  region  of  nonzero  vorticity.  In  addition,  all  boundary  conditions,  and  in  par¬ 
ticular  the  farfield  boundary  condition,  are  naturally  incorporated  in  the  formula¬ 
tion. 
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APPENDIX  A 


DERIVATION  OF  VECTOR  INTEGRAL  IDENTITY 


One  may  start  by  noting  that  the  vector  identities 

Vx(Gxu)  =  (u-V)G-(G-V)u  +  G(V-u)-u(V*G), 
V(G-u)  =  (G  *  V)u  +  (u  •  V)G  +  G  X  (V  X  u)  +  u  x  (V  x  G) 


may  be  added  with  the  result  that 

Vx(Gxu)  +  V(G-u)  =  2(u  •  V)G  +  G(V  •  u)  -  u(V  •  G) 

+  G  x  (V  x  u)  +  u  x  (V  x  G).  (A-2) 


If  one  integrates  this  expression  over  the  volume  of  the  domain  and  assumes  that 
the  functions  under  examination  are  sufficiently  well  behaved  for  the  integrals  to 
exist,  then  one  may  apply  the  curl  and  gradient  theorems, 


rrr  rr 

V xFdV  =  <D>  nxFdS, 


jj 


VfdV  =  (H)  fndS, 


(A-3) 


to  the  volume  integral  of  the  left-hand  side  of  the  expression  to  find  that 

[V  x  (G  x  u)  +  V(G  •  n)]dV  =  <jj>  [n  x  (G  x  u)  +  (G  •  u)n]dS.  (A-4) 


JJJ 

V 


JJ 

s 


However,  it  can  readily  be  shown  that 


nx(Gxu)  =  (n  •  u)G  -  (n  •  G)u, 
(G  ■  u)n  =  G  x  (n  X  u)  +  (n  ■  G)u, 


(A-5) 


A-l 


so  that 


[Vx(Gxu)  +  V(G  •  u)]dV  =  #[(  n  •  u)G  +  G  x  (n  x  u)]  dS.  (A-6) 


One  is  then  left  with  the  expression 


{(n  •  u)G  +  G  x  (n  x  u)}  dS  =  Jjj  [2(u  •  V)G  +  G(V  •  u)-  u(V  •  G) 

*v 

+  G  x  (V  x  u)  +  u  x  (V  x  G)]dV.  (A-7) 


The  divergence  theorem  may  then  be  used  to  show  that 


(u  •  V)GdV  =  (ip  (n  •  u)GdS  -  (V-u)GdV, 


(A-8) 


so  that  equation  (A-7)  becomes 


[(n  •  u)G  -  G  x  (n  x  u)]dS  =  {G(V  •u)+u(V*  G) 


■  G  x  (V  x  u)  -  u  x  (V  x  G)}  dV.  (A-9) 


Now  if  one  lets 


G  =  — i —  +  VH, 
r(r  +  e)2 


(A-10) 


where  H  is  some  function  that  is  regular  in  the  fluid  domain,  then,  in  the  limit  as  e 
approaches  zero,  one  finds  that 


u(V  •  G)dV  =  -  4'7ru, 


u  x  (V  x  G)dV  =  0. 


(A-ll) 


A-2 


Hence,  equation  (A-9)  may  be  written  as 


4-7711  =  -  #[(  n  •  u)G  -  G  x  (n  x  u)]dS 


+  [(V  •  u)G  -  G  x  (V  x  u)]dV. 


(A- 12) 


In  general,  if  the  field  point  is  taken  to  be  in  the  domain,  on  the  boundary  of 
the  domain,  or  in  the  complement  of  the  domain,  then  the  expression  becomes 


4ir  in  D 

277  on  S  ►  u  = 
0  in  Dc  „ 


[(n  •  u)G  -  G  x  (n  x  u)]dS 


[G  x  (V  x  u)  -  (V  •  u)G]dV. 


(A- 13) 


If  G  is  taken  to  be  the  free-space  Green’s  function 


(A-14) 


then  equations  (A- 12)  and  (A-13)  become 


4l7U  =  - 


(n-tt)r  .  rx(nxu)ldg 

r3  r3  J 


rX(V  xu)  (V-u)r 
r3  "  “r3 


]dV, 


(A-15) 


2ttu  =  -  ©  ffou)r  -  rX<nXu)ldS 

Jj  L  T3  T3  J 


rX(Vxu)  _  (V 


I2*l< 

r3  J 


(A- 16) 
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APPENDIX  B 


DERIVATION  OF  SCALAR  INTEGRAL  IDENTITY 


One  may  begin  this  derivation  with  Green’s  second  identity  (see  reference  3): 
ri.n2A.\A _  <il\  ( u.  rL 


($V2G  -  GV20)dV  =  #  (<t>  -  G  ^  )dS. 


(B-l) 


Now  let 


G  =  ,— ±_  +  H, 

(FT7) 


(B-2) 


where  H  is  some  function  that  is  regular  in  the  fluid  domain.  Then  it  can  be  shown 
that 

“So  ffl  (B-3) 


lr(rh> 

V 


so  that,  in  the  limit  as  e  approaches  zero,  one  finds  that 

rrr 

V20G  dV, 


JJ 

s 


(B-4) 


and,  in  general,  if  the  field  point  is  in  the  domain,  on  the  boundary  of  the  domain, 
or  in  the  complement  of  the  domain,  this  becomes 


4 tt  in  D 
2rr  on  S 
0  in  Dc 


1 

7*  /  ,  _ . 

'ti 

II 

>(lG-*i)ds- 

i 

V20GdV. 


(B-5) 


v 
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If  G  is  taken  to  be  the  free-space  Green’s  function 


G  r’ 


(B-6) 


then  equations  (B-4)  and  (B-5)  take  on  the  usual  forms  of  Green’s  third  identity: 


4^=#fr4(;)]dS-  dV’ 


(B-7) 


2^  =  |[lr4(7)]dS- 

JJ  JJj 


(B-8) 
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APPENDIX  C 


PROOF  OF  WALL  VORTICITY  VECTOR’S  TANGENCY 


For  any  two  orthogonal  unit  vectors  tangent  to  the  boundary  t;  and  si}  the  unit 
vector  normal  to  the  boundary  may  be  expressed  as 

ni  =  eijktjsk  5  (C-l) 


where  £jjk  is  the  alternating  tensor.  Since  the  vorticity  is  defined  as  (see  Jeffreys 
(reference  5)) 


(C-2) 


the  normal  component  of  the  vorticity  at  the  wall  may  be  written  as 


(C-3) 


Each  of  the  last  two  terms  represents  a  directional  derivative  tangent  to  the 
boundary.  If  one  assumes  that  the  boundary  motion  is  rigid,  then  on  the  boundary 

Uj  =  Uj  +  £ijkXjflk  ,  (C-4) 

where  Uj  is  the  translational  motion  of  the  boundary  and  ft;  is  its  rotational  motion 
about  the  origin  of  coordinates.  With  this  representation  one  finds  that 


C-l 


3uj 


J3xj 


—  ejjktjftk, 

=  ni(skfik)  -  Sj(nkflk) , 


(C-5) 


and 


Sjgx.  £ijksj^\  > 

'  =ti(nkak)-ni(tknk),  (C-6) 


so  that 


n^COj  Sk(tj£kjmflm)  ■  tj(Sk£jkmIlm): 

^£kjm^k^j^,m , 

=  -2nmOm.  (C-7) 


Thus,  for  rigid  boundary  motion,  the  normal  component  of  the  vorticity  may  be 
determined  from  the  velocity  boundary  conditions. 


The  cross  product  of  the  normal  and  the  vorticity  at  the  wall  may  be 
determined  in  a  similar  manner.  The  cross  product  is  defined  as 


£ijknjt0k 


(C-8) 


If  one  employs  equation  (C-2),  this  expression  becomes 


_  3um 

£ijknjWk  —  ^ijk^klm^j  > 


3U„ 


\dxi  3xj  / 


(C-9) 


C-2 


so  that 


9U; 


3ui 


3u: 

- - is,- 

k3xk 

3ui  \ 
f  nv. — -n;  1 

3uj 

'  nj^’ 

k3xk  7 

.  3v 

>k9xk 

v  /  3uj 

3uk 

nkIlj  3xj 

nj[(n_j  skflk  -  Sj  nknk)ti  +  (tjnkak  -  njtkflk)si] 

/  0U;  3uk  \ 

■(n^-n^Di)’ 

(9ui  9uk  \ 
nj  -gp  -  nknj  -^:nij ; 


hence,  the  components  of  the  cross  product  are 

3uj 

ti(eijknj(ok)  =  skOk  -  tjnj  — , 

.  .  _  3ui 


(C-10) 


(C-ll) 


Since,  for  rigid  boundary  motion,  the  cross  product  of  the  boundary  normal  vector 
and  the  boundary  vorticity  contains  derivatives  of  the  velocity  in  the  direction 
normal  to  the  boundary,  this  cross  product  cannot  be  determined  from  the 
boundary  conditions  alone. 
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